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ABSTRACT 

A number of recent results from X-ray observations of Active Galactic Nuclei 
involving the Fe line (reduction of line variability compared to the X-ray continuum 
variability, the X-ray "Baldwin effect") were attributed to a presence of a hot, ionized 
skin of an accretion disc, suppressing emission of the line. The ionized skin appears as 
a result of the thermal instability of X-ray irradiated plasma. We test this hypothesis 
by computing the Thomson thickness of the hot skin on top of the aPtot Shakura- 
Sunyaev disc, by simultaneously solving the vertical structure of both the hot skin and 
the disc. We then compute a number of relations between observable quantities, e.g. 
the hard X-ray flux, amplitude of the observed reprocessed component, relativistic 
smearing of the line, the r.m.s. variability of the hard X-rays. These relations can 
be compared to present and future observations. We point out that this mechanism is 
unlikely to explain the behaviour of the X-ray source in MCG-6-30-15, where there is 
a number of arguments against the existence of a thick hot skin, but it can work for 
some other Seyfert 1 galaxies. 



Key words: accretion, accretion disc 
Seyfert - X-ray: galaxies 



instabilities - galaxies: active - galaxies: 



1 INTRODUCTION 

Spectral features originating as a result of reprocessing of 
energetic X-ray radiation by cool, optically thick plasma 
around accreting compact objects can be extremely help- 
ful in constraining the geometry of accretion flows. Indeed, 
one of the strongest evidences for the existence of accretion 
discs very close to central black holes in active galactic nuclei 
(AGN) came from observations of broad profiles of the fluo- 
rescent iron Kc, line near 6.4 keV in many Seyfert 1 galaxies, 
the ASCA result from MCG-6-30-15 (Tanaka et al. 1995) be- 
ing the best known case. The broad line profile is unlikely to 
be explained as a result of Comptonization in transmission 
through a plasma cloud (Fabian et al. 1995), so broadening 
as a result of Doppler effects and gravitational redshift seems 
to be the main mechanism (although Comptonization in re- 
flection can contribute to the broadening; Karas et al. 2000) . 
Quantitative estimates of the distances of the cold plasma 
from central black holes based on the line profiles indicate an 
intrinsic dispersion of parameters among the sources (Nan- 
dra et al. 1997; Done, Madejski & Zycki 2000), but generally 
yield the distance < lOOrg (rg = GM/c^). 
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Data analyzed in most cases are integrated over rel- 
atively long periods of time (~days) in order to decrease 
statistical errors. They are therefore averages over many 
dynamical timescales, tdyn ~ , and correspond to 
time/space average geometrical/physical conditions in the 
region where the line is produced. Studies of time variability 
of the line can add a new dimension to inferences based on 
averaged spectra, as they can potentially resolve the struc- 
ture of the X-rays production region. Reverberation analysis 
(investigations of the response of the line to variability of the 
continuum driving its emission; Blandford & McKee 1982) 
can provide independent constraints on the geometry of the 
emission region (e.g. Reynolds et al. 1999). 

Indeed, a number of surprising results were obtained in 
those cases where a good quality of data enabled studies of 
the line variability. Generally, the line absolute flux appears 
to be less variable than the flux of the driving continuum. 
The best example here is the long campaign of Rossi X- 
ray Timing Explorer (RXTE) observations of MCG-6-30- 
15 analyzed in Reynolds (2000) and Lee et al. (2000). The 
r.m.s. variability in the line band was reduced to ~ 0.8 x 
that of the continuum, due to a combination of constant line 
flux and flux-correlated spectral changes. Similarly, Done et 
al. (2000) analyzed RXTE observations of another Seyfert 
1 galaxy IC 4329a and found that the absolute amplitude 
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of the reprocessed component stayed roughly constant as 
the source brightened, i.e. the relative amplitude decreased 
somewhat. Chiang et al. (2000) found constant line flux in 
NGC 5548, even though the 2-10 keV continuum flux varied 
by a factor of 2. 

However, the lack of response of the line integrated flux 
to the continuum variability does not mean that there are 
no changes at all. The profile of the line clearly responds to 
the continuum variability. Again, the most striking example 
is MCG-6-30-15, where the line profile showed quite dra- 
matic changes (Iwasawa et al. 1996) : it was narrow during a 
high flux state but it was very broad when the source's flux 
decreased. Interestingly, the total flux of the line photons 
seems to have been constant in time (see fig. 6 in Iwasawa 
et al. 1996). 

A suggestion that the line flux can be constant due 
to ionization of the reprocessor varying in relation to the 
primary flux were made by Reynolds (2000). However, no 
physical model was proposed for the flux-ionization de- 
pendence. Nayakshin, Kazanas & Kallman (2000; hereafter 
NKK) pointed out that a reduction of the Fe line vari- 
ability can be expected if the disc develops a hot, ionized 
layer where the line production is inefficient. Development of 
such an ionized layer is a necessary consequence of the ther- 
mal instability operating in X-ray irradiated plasma (Krolik, 
McKee & Tarter 1981; Rozahska & Czerny 1996, hereafter 
RC96) . Since the thickness of the hot layer and its radial ex- 
tent are correlated with the illuminating flux, the efficiency 
of the Kc line production can indeed be anti-correlated with 
the source X-ray luminosity. Moreover, changes to the ra- 
dial location of the line emission region will influence the 
line profile. 

The same idea has been proposed by Nayakshin 
(2000a,b) to explain the X-ray "Baldwin effect" (Baldwin 
1977; Iwasawa & Taniguchi 1993; Nandra et al. 1997), i.e. 
smaller equivalent width of the Kq, line for X-ray brighter 
sources. The two effects (the Baldwin effect and the lack of 
reverberation signatures) can be simply considered as the 
same effect looked at statistically in a sample of objects and 
in a time history of a single object, respectively ("ergodic 
theorem"), although the two phenomena may differ in de- 
tails: X-ray variability of a single source is probably related 
to variable output of some magnetic activity, while differ- 
ent (averaged) luminosities of sources are probably due to 
different mass accretion rates. 

In this paper we consider and test this model further, 
based on computations of vertical structure of illuminated 
accretion discs in hydrostatic equilibrium (RC96; Rozahska 
1999, hereafter R99; see also NKK). The plan of the paper is 
as follows: in Section |^ we sketch the physics of the thermal 
instability responsible for the two-layer structure of irradi- 
ated accretion discs, in Section ^ we present our method of 
computing the thickness of the hot layer and results of such 
computations. Section ^ contains discussion of consequences 
of the presence of the hot layer for variability properties of 
accreting objects, while Section |5| contains similar discussion 
on spectral properties. The results are discussed in Section pi 



2 PHYSICAL PICTURE 



2.1 Thermal instability of irradiated plasma 

Our physical picture is based on vertical structure of illu- 
minated accretion discs. As a result of illumination the top- 
most disc layer can be strongly ionized, but the level of ion- 
ization decreases with depth into the disc. In calculations 
performed with the assumption of constant density (Ross 
& Fabian 1993; Zycki et al. 1994), the ionization level at 
the disc surface is given by ^ = Fx/ric {Fx is the ffux of 
illuminating X-rays and Uc is the electron number density) 
and it smoothly decreases inwards. Adjusting any of these 
parameters it is possible to obtain any specified ionization 
stage. The picture is rather different when the condition of 
hydrostatic equilibrium (or in fact any condition on gas pres- 
sure rather than its density) is additionally imposed in the 
illuminated disc (RC96; NKK). Solutions corresponding to 
certain range of temperatures (and/or ionization parameter 
H = Fx/{cPgas), more appropriate than ^ for characterizing 
solutions with pressure constraints) are thermally unstable 
(Field 1965; Krolik et al. 1981). The condition for instability 
(Field 1965), 

where £ is the heat-loss function, can be expressed as 

(§)<»■ 

The unstable solutions are thus those located on the part 
of the T-H diagram with a negative slope (Krolik et al. 
1981). The physical mechanism for the instability is the 
weak dependence of bremsstrahlung cooling on tempera- 
ture, jbroms oc T^^^, which cannot balance the photoion- 
ization heating on heavy elements caused by hard X-rays. 
On the lower stable branch at T ~ 10* K the plasma is 
cooled by atomic line emission, while on the upper branch, 
at T ~ 10^-10* K, Compton heating and cooling are the 
dominant processes. 



2.2 Vertical structure of irradiated discs 

In an illuminated accretion disc in hydrostatic equilibrium 
a movement along the T-H curve corresponds to chang- 
ing the depth into the disc. As a result of the instability, 
a broad range of intermediate ionization of iron is inacces- 
sible to a stable solution. The uniqueness of the solution 
is assured only when the basic equations: energy equation, 
hydrostatic equilibrium and radiative transfer equations are 
supplemented by the heat conduction equation (Maciolek- 
Niedzwiecki, Krohk & Zdziarski 1997; R99). The transition 
from the hot to cold solution branch is very sharp and 
its location can be obtained from the general condition of 
radiative-conductive equilibrium (Rozahska 2000) . The con- 
dition gives the temperature at the transition point 

7 

^trans — Tj^-^s, (3) 

where Ts is the surface temperature. For purely radiation 
heated disc atmosphere Tg = Tic, where Tic is the inverse- 
Compton temperature (McKee & Begelman 1990; Rozahska 
& Czerny 2000a). When the aP heating is also included, Ts 
is somewhat higher than Tic (Rozahska 2000). In a good 
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approximation, the transition occurs where the upper, sta- 
ble branch of the H-T curve changes to the middle, un- 
stable branch (RC96, R99, NKK), The resuhing structure 
is then basically two-layered: the upper, hot layer (HL) and 
the lower, cold layer, although there is a thin transition layer 
in between the two (R99). For a sufficiently hard illuminat- 
ing spectra (i.e. high Tic) iron can be considered completely 
ionized in the hot layer, while it recombines to below FeXI 
in the cold layer (NKK; see also fig. 2 in Zycki & Czerny 
1994). 



3 THE THICKNESS OF THE HOT LAYER 

The model parameter crucial for spectral/timing predictions 
is the thickness of the hot layer, rhot- The larger the thick- 
ness, the smaller the fraction of primary radiation penetrat- 
ing to the cold layer and giving rise to the usual reprocessed 
component. The reprocessed component is further Comp- 
tonized as the photons escape through the hot layer. The 
result (at least in the limited energy band corresponding to 
e.g. Ginga or RXTE data) is a reprocessed component with 
amplitude reduced compared to the usual geometrical fac- 
tor n/47r (Eq. where Q. is the solid angle subtended by 
the reprocessor from the X-ray source (Done & Nayakshin 
2001). The amplitude of the iron spectral features is further 
reduced by the Comptonization. 

Two most important parameters determining the thick- 
ness of the hot layer are: the strength of irradiation com- 
pared to internal disc emission, Fx/Fu, and the ratio of 
gravity at the base of the layer to X-ray radiation pressure 
(RC96, NKK). Following NKK we parameterize the latter 

by 

where 

/x = :t—, (5) 

moC^ 

where = p/rm ~ 1.9 for cosmic abundance. Obviously 
then, the structure of the hot layer has to be solved simulta- 
neously with the structure of the cold disc, since A depends 
on the disc thickness, Hd- NKK used vertically averaged so- 
lutions of Shakura & Sunyaev (1973; hereafter SS) for the 
disc in their computations of rhot. They assumed that to- 
tal pressure, Pgas -f -Prad, should be continuous across the 
boundary between the two layers, and obtained both Pgas 
and Prad individually continuous across the boundary. 

In this paper we explicitly solve equations of vertical 
structure of the cold disc. We do this for the aPtot SS disc, 
i.e. a disc with energy generation proportional to total pres- 
sure. We define the input parameters of the computations 
in such a way as to make the connection between accretion 
discs models and observable quantities obvious. First, we de- 
scribe calculations at a given radius and then discuss radial 
dependence of rhot. 



3.1 Vertical structure at a given radius 

The structure of the HL at a given radius is computed by 
combining the method of RC96 and R99 with computa- 
tions of the vertical structure of X-ray illuminated discs by 



Rozanska et al. (1999). For the HL we follow closely the 
method of RC96 and R99, with the important simplifica- 
tion of neglecting thermal conduction. The reason is purely 
practical, as efforts to fully combine proper photo-ionization 
computations of the hot layer with vertical disc structure 
have only just begun (Dumont, Abrassart & Collin 2000; 
Rozanska et al., in preparation). This leaves us with an im- 
portant problem of selecting proper solution in the zone of 
instability, where the system of equations has more than one 
solution for temperature. We adopt a simple prescription 
and select the highest value of T as the solution. Comparing 
solutions with and without conduction in RC96 and R99 one 
sees that our procedure may overestimate somewhat the to- 
tal thickness of the hot and transition layers (i.e. the depth of 
the point where T = Tctt). However, the thickness of the hot 
layer alone is not afi'ected by our neglecting the thermal con- 
duction. The spectrum of illuminating radiation is assumed 
to be a power law with a cutoff, Fe oc E~^'^^e~^^^'' , param- 
eterized by the photon spectral index F and cutoff energy. 

Equations of vertical structure of the cold disc are the 
same as in Rozanska et al. (1999). We assume that a certain 
fraction of gravitational energy, ^, is dissipated within the 
disc (but the disc transports all the angular momentum, see 
e.g. Witt, Czerny & Zycki 1997). The remaining fraction, 
1 — ^, is dissipated in an active corona and converted to 
hard X-ray radiation, of which a fraction O.Stj illuminates 
the disk, i.e. 



Fx = 0.5r?(l -O^lot, 
where 

3 GMM ( 
Ptot = G(r) = -^:^ 1 




(6) 



(7) 



is the gravitational energy dissipation per unit area of the 
disc surface. Value of 77 = 1 would thus correspond to a sta- 
tionary illumination by an X-ray producing corona above 
the disk, assuming isotropic emission (see e.g. Haardt & 
Maraschi 1991). We will assume > 1, corresponding to 
a possible enhancement of the X-ray flux. Such an enhance- 
ment can be expected if the X-rays are produced by mag- 
netic flares, due to 1) the "charge" time of a flare being 
longer than the "discharge" time, 2) the collecting area for 
the flare energy being larger than the illuminated area (e.g. 
Haardt, Maraschi & Ghisellini 1994). 

We solve the structure of the HL assuming a certain 
initial geometrical thickness of the total system, //tot- Four 
parameters, computed at the bottom of the HL, are then 
used to modify the equations of the disc structure and their 
boundary conditions: 

• Location of the bottom of the HL (height above the 
equatorial plane) gives the initial thickness of the cold disc, 

• The hard X-ray flux attenuated by the HL contributes 
to heating of the disc. 

• The temperature at the top of the cold disc is com- 
puted from the Eddington approximation, T — Toff [1/2 -|- 
(3/4) r]^''*, taking r to be equal to the optical depth of the 
HL, and Toff from both the internal energy dissipation and 
the irradiation, 

• Gas pressure at the bottom of the HL determines the 
gas density at the top of the disk. 
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Figure 1. Vertical structure of tlie liot layer on top of an accretion disc at r = lOr-g, for M = 10 Mq, m = 0.05, a = 0.1, 5 = 0.5 
and rj = I. Plotted are the gas temperature, gas pressure and density. Sharp transition between the upper, hot, ionized layer and cool 
disc is observed. The thickness of the hot layer (from Tg Kb Tjc to 7/15 Ts) is rjjot ^ 0.1. The non-monotonicity of Pgas and p is due to 
contribution to energy transfer from vertical convection, and contribution from radiation pressure. 



Iterations of the vertical structure of the disc now fol- 
low, to compute the disc thickness for given disc dissipation 
^Ftot, with the boundary conditions as above. Usually, the 
converged thickness of the disc is different than the initially 
assumed one. Therefore the final solution (at a given radius 
r and for an assumed central mass M, mass accretion rate 
m = M/A/Edd, and viscosity parameter, a) is obtained by 
iterating i/tot until a consistent solution is found. As already 
mentioned above, we define rhot as the thickness of the layer 
where temperature changes from Ts down to (7/15) Ts. 

As a reference, we compute Thot for parameters appro- 
priate for a Seyfert 1 galaxy: M = 10^ M©, m — 0.05, 
r = lOrg (close to the peak in energy dissipation), a — 0.1, 
^ = 0.5 and = 1 (i.e. no enhancement), F — 1.9 and 
Ec = 150 keV. 

Figure |l| shows an example of the vertical structure of 
the disc with the hot layer. The Thomson thickness of HL is 
rather small, rhot ~ 0.1, irrespectively of a. The "gravity pa- 
rameter" is A 7 in this case. Our value of Thot is thus very 
similar to what NKK obtain for this value of A (see fig. 4 in 
NKK), demonstrating that, despite some approximations in 
our procedure, our results are accurate. 

The main parameter, A, is a non-unique function of 
other parameters of the system. For given m and r, A is 
determined by a combination of ^ and i]: both parameters 
influence the disc thickness and the irradiation flux. We have 
therefore extended our reference computations to obtain rhot 
on a grid of £, and rj. Results are plotted in Figure |^, where 
we show a number of quantities characterizing the disc or 
the spectrum. 

The amplitude of the observed "cold" reflection 
(energy-integrated albedo) is estimated using 



R 



btr(Thot)]^ 
1 +Prf(Thot) 



(8) 



where ptr(rhot) and Prf(rhot) are hard X-rays transmission 
and reflection probabilities for the HL thickness Thot- They 
are defined as ratios of the transmitted and reflected flux 
to the incident flux, respectively, for a purely scattering at- 
mosphere. The factor p^r appears in the numerator of equa- 



tion (|8|), because the reflected photons have to go through 
the HL twice. Additionally, a fraction prf of the incident pho- 
tons is elastically back-scattered off the HL, and, by our as- 
sumption, simply contributes to the primary continuum pho- 
tons initially directed towards an observer. In the limiting 
case of no HL (rhot 0), ptr(-rhot) 1, while Prf(Thot) 
and i? — > 1, i.e. 7? represents then the solid angle of the 
reprocessor, normalized to 2tt. We compute ptr(rhot) and 
Prf(''"hot) using a Monte Carlo simulations of Compton scat- 
tering. This is a rather approximate estimate of the angle 
averaged amplitude, since it assumes no photo-absorption 
in the HL, i.e. no contribution to the reprocessed spectrum 
from highly (but not completely) ionized elements. It thus 
assumes a rather hard illuminating spectrum (F < 2) giving 
high Tic. A more accurate procedure of estimating R was 
employed by Done & Nayakshin (2001), who fltted spectra 
computed by NKK by a "cold" reflection model in the con- 
text of Ginga data. For softer illuminating spectra, F <; 2, 
the reprocessed component contains strong spectral features 
due to H- and He- like Fe ions and it cannot be approximated 
by "cold" reprocessing (NKK; Nayakshin 2000b). 

The ratio -Fhard/^soft is the observed ratio of fluxes, i.e. 
the hard X-ray flux 



fhard = -Fx(l +Prf) 



(9) 



contains a contribution from photons backscattered in the 
HL towards an observer. The soft flux, -Fsoft, contains a con- 
tribution from the disc internal emission, _Fd = C^tot, and 
hard photons thermalized in the cold disc. 



Tsoft =Fo + {l- a)ptr-Fx, 



(10) 



with a ft; 0.2 being the energy-integrated X-ray albedo of 
cold matter. We note that the above relations of -Fhard and 
-Fsoft with Fx assume implicitly that there are no significant 
changes of the source geometry (e.g. the height of the flare 
above the disc) , so that the variability of _fx is solely due to 
the variable output of the magnetic activity. 

The structure of the disc has to be rather different from 
the standard SS solution, if the hot layer is to be Thomson 
thick(ish), Thot ^ 1- For half of the total energy dissipated 
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Figure 2. Results of computations of the structure of tlie liot layer on top of an accretion disc at r = 10 rg, for M = 10^ Mq , m = 0.05 
and viscosity parameter o = 0.1. The panels show contours of: i^x/^D ^ the ratio of the illuminating X-ray flux to the disc flux due to 
internal dissipation; A - the gravity parameter (see equation tjjqj - the Thomson thickness of the hot layer; the effective amplitude 
of the cold, reprocessed component (equation B); i^hard/^soft ~ the ratio of the observed hard X-ray flux to the soft, thermal flux. 



in the active regions, ^ = 0.5, the illumination flux has to 
be further enhanced by a factor of 77 « 40 relative to the 
stationary situation in order to obtain rhot = 1- This cor- 
responds to the ratio of illumination flux to disc internal 
flux, Jx/Jb = 0.5?7(1 - 0/C ~ 20. For larger fraction of 
energy dissipated in the magnetic corona, i.e. smaller ^, the 
required ratio -Fx/Jb is even larger (e.g. Fk/Fd ~ 60 for 
5 = 0.1). This appears to be due to the disc thickness, H^, 
increasing with ^ more slowly than linearly, and thus for 
A oc (Fo/ Fx){Hd/(,) to remain constant, Fx/Fo has to in- 
crease when ^ decreases. 

On the other hand, the presence of the HL will be appar- 
ent in observations even if its thickness is rather smaller than 
Thot = 1. According to our approximate formula (equation^ 
reflection amplitude gets reduced to 0.8 for rhot ~ 0.25. 
Quality of current data is usually sufficient to distinguish 
R = 0.8 from _R = 1 in typical X-ray spectra of AGN (e.g. 
Done et al. 2000). 




10 100 

Radius (R ) 



Figure 3. Thomson thickness of the hot layer as a function of 
radius for r] = 3 and three values of as labelled. Other param- 
eters: M = lO'^ Mq, m = 0.05 and a = 0.1. 



3.2 Radial structure 

The Thomson thickness of the HL changes with radius, even 
for a constant Fx/Fo, since A is not constant with r. We 
computed radial dependencies of rhot for rj — 3 and three 
values of ^ = 0.01,0.3 and 0.6. Values of other parameters 



were as before: M = 10^ M0, m = 0.05, a — 0.1. Results 
are plotted in Fig. ^. 

The thickness is largest at r « 10 rg. It then drops 
slowly, roughly in accord with the formula of Nayakshin 
(2000a). A consequence of that particular rhot('') distri- 
bution is a relative suppression of highly relativistically 
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Figure 4. Upper panel: contours of observed amplitude of the 
reprocessed component (solid curves, labelled), and contours of 
observed hard X— ray flux (dashed curves). Each contour of Fhard 
means an increase of Fhard by a factor of 2. Lower panel: pro- 
jected dependence of i? vs. i^hard- Clearly, the dependence is much 
weaker than linear, i.e. it is not possible for the absolute normal- 
ization of the reproce^ed component, N = R X Fhard to remain 
constant. See Figured for parameters of computations. 



smeared component of the reprocessed spectrum. On the 
other hand, the reprocessed photons that do manage to dif- 
fuse from below the hot layer close to the central BH, suffer 
Comptonization smearing, which in practice may be difficult 
to distinguish from relativistic smearing (see Sec. [5.2| ). 

Obviously, in the magnetic flare scenario there is no 
reason for changes of ^ and rj to be perfectly correlated on 
all radii. One can thus imagine rather different values at 
different radii, leading to Thot('") rather different than those 
plotted in Fig. |^. On timescales long compared to tdyn one 
could perhaps obtain time average ^(r) and ri{r) employing a 
particular variability model, e.g. that of Poutanen & Fabian 
(1999). 



3.3 Dependence of Thot on X— ray flux 

The qualitative relation between the thickness of the HL 
and the illuminating X-ray flux is obvious: the stronger the 
X-ray flux the thicker the HL. In this Section we present 
the quantitative relation from our model. 

In Figure ^ we again plot contours of amplitude of re- 
flection R in the ^-ry plane, together with contours of the 



observed X-ray flux, Fhard, as well as the projected depen- 
dence -R--Fhard- For 7? close to 1, doubling Fx results in only 
modest decrease of R, by 10-20 per cent. For smaller R the 
process is more effective. But still, to go from e.g. R — 0.8 
to R — 0.4 requires an increase of -Fhard by a factor of 8 
(when the movement is along the ^ axis, for « 4). Going 
from R = 0.6 to R = 0.3 requires that Fhard increases at 
least 4 times, again approximately along the ^ axis. In no 
case does it seems possible to obtain an inverse proportion- 
ality of Fhard and R, so that R x Fhard ~ const. Thus, these 
results seem to rule out the possibility that the absolute flux 
of the Fe Kq, line may remain constant despite continuum 
variability, although certain reduction of the line variabil- 
ity amplitude (compared to that of the continuum) can be 
expected (see next Section). 

We note that a change of Fx is more effective in chang- 
ing Thot , if Fx increases or decreases due to changing fraction 
of energy dissipated in the disc, ^, rather than rj. The reason 
is that e.g. removing the energy dissipation from the disc to 
the active corona (i.e. decreasing ^) not only increases the 
X-ray illuminating flux, but it also decreases the disc thick- 
ness, thus decreasing the gravity. 



4 REDUCTION OF VARIABILITY 
AMPLITUDE OF THE FE LINE 

4.1 Method 

Based on results of previous Section we expect the variability 
of the Fe Kq, line flux to be reduced compared to the ampli- 
tude of variability of the driving continuum, even though it 
does not seem possible to obtain absolutely constant flux of 
the line. Here we estimate the reduction of the r.m.s. vari- 
ability of the Fe line, compared with an assumed r.m.s. 
variability of its driving continuum. To this end we will use 
results of our calculations of rhot(^,??) at r = lOr-g. In this 
way we will obtain the upper limit to the considered effect, 
since Thot(f') reaches maximum at r ~ 10 rg, for a given pair 
(C,r?) (see Fig. |). 

First, we construct a simulated light curve of hard X- 
ray continuum by summing Fourier components with ran- 
dom phases (e.g. Tsonis 1992), 



C{U) = J2 VpUT) cos [27rf,U + m 



(11) 



for an adopted form of power spectral density (PSD), P(f) 
For the PSD we adopt 

1 



[1 + (///c 



(12) 



This provides a good description of hard X-ray PSD of e.g. 
MCG-6-30-15, with the cutoff frequency /c = 10"^ Hz, the 
slope p = 1.4 and the normalization corresponding to r.m.s. 
variability of 20 per cent (Green, McHardy & Lehto 1993; 
Czerny & Lehto 1997; Yaqoob et al. 1997; Nowak & Chiang 
2000; Reynolds 2000). We take these parameters to be rep- 
resentative for Seyfert 1 galaxies. The light curve is scaled to 
have a given mean value C. Next, the light curve in energy 
band containing the line photons is constructed accord- 
ing to 



L{U) = (1 - Ai)C{U) + AiCiU)R[C{U)], 



(13) 
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Figure 5. Examples of light curves of the observed hard X-ray 
continuum (open circles) and the Fe Kr> l ine (solid squares), sim- 
ulated from the PSD defined in Sec The assumed parame- 
ters gave the effective amplitude of the reprocessed component 
R{C) ~ 0.75. In effect the line flux is suppressed at almost all 
times and reduction of the line r.m.s. variability is significant, at 
6 per cent (see Fig. 



where Ai is the fraction of line photons contributing to to- 
tal counts in the considered energy band, while R[C{t)] is 
the relative amplitude of the reprocessed component as a 
function of the illuminating hard X-ray flux. Writing equa- 
tion (|l|) we have assumed no time lag in the response of the 
Fe Ka line emission after the continuum variations, consis- 
tent with observations (e.g. Reynolds 2000). The amplitude 
R[C{t)] is computed as follows: we choose a value for the ef- 
fective amplitude of the reprocessed component at the mean 
flux level, R{C). We then assume that R{C) corresponds to a 
pair (^0 , Jyo) , which is located on a trajectory in the £,-t] plane 
that the system travels while changing its luminosity. The 
trajectory describes how changes of ^ and rj contribute to 
the variability. Since our understanding of variability mech- 
anisms is rather limited, we will simply test a number of 
arbitrary trajectories, to get an idea of available range of re- 
sults. As discussed in previous Section, variations along the 
^-axis give largest possible change of R for a given change 
of -Fx, therefore our chosen trajectories have a significant 
component along the ^ axis. 

Next, we construct the light curve of the observed 
hard X-ray continuum, which has a contribution from pho- 
tons backscattered in the hot layer, i.e. F{t) — C{t){l + 
Pii{Thot[C{t)]}}. Finally, the r.m.s. variability is computed 
for F{t) and L{t) according to 



0.4 



R(C) 



.(C) = = 



\ 



(14) 



4.2 Results 

First we set Ai = 1 i.e. we consider the light curve of the Fe 
Kq line alone, without any continuum photons. An exam- 
ple of the simulated light curves is plotted in Fig. H. They 



Figure 6. Amplitude of variability of the Fe Kc« line flux is re- 
duced compared the that of the driving continuum, due to ap- 
pearance of a hot, ionized layer with thickness increasing with 
increasing illuminating flux. This is plotted here as a function 
of the assumed effective amplitude of reflection at the mean flux 
level. Upper panel shows reduction of the r.m.s. variability as- 
suming that there are no primary continuum photons in the line 
energy band. Lower panel shows analogous reduction assuming 
that the line photons contribute 15 per cent to the considered en- 
ergy band. Different symbols denote different trajectories in the 
£^—r] plane traveled by the varying system. The r.m.s. va riab ility 
of the driving continuum was 20 per cent. See text (Sec 4.2) for 
further explanations. 



were obtained assuming that the mean flux of the intrinsic 
continuum, C, corresponds to fp = 0.5 and -qo — 1 (which 
yields R{C) « 0.75, see Fig. W), and that the variability 
is only due to changing ^. The line flux is reduced at all 
times, since R{C) is significantly lower than 1, and its r.m.s. 
variability is reduced by « 5 per cent compared with the 
r.m.s. variability of its driving continuum, C{t). The ob- 
served hard X-ray continuum, F{t), contains a contribution 
from photons backscattered in the hot layer, which increase 
its r.m.s. variability by ~ 1 per cent. Thus, the difference of 
r.m.s. variabilities between the line and the hard continuum, 
f7"var(-F) — f'"var(i) is in this example ~ 6 per cent. 

In real spectra the contribution of line photons to the 
5-7 keV band is only Ai ~ 0.15, for R — 1. The reduction 
of r.m.s. variability, Auvar = crvar(-F) — o"var(£'), is therefore 
much weaker than in previous case, and in the above exam- 
ple, only ~ 2 per cent. 

The difference Aavar is plotted in Fig. |^ as a function of 
R{C), for both j4i = 1 and Ai = 0.15. The general trend is 
that the larger the R{C) is, the smaller the avar(-F') — o"var(i) 
is. Consider for example R{C) « 1, i.e. at the mean ffux 
level the irradiation flux is relatively weak and in conse- 
quence the hot layer is rather thin, rhot < 0.1. Then R is 
reduced only when the instantaneous flux exceeds the mean 
flux, but R saturates at 1 whenever C'{t) < C. This then 
gives smaller reduction of r.m.s. variability compared to a 
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Figure 7. Panel a): the amplitude of reflection — spectral index correlation predicted by the model for a number of fis values (equation |l6[). 
Only models with ^ < 0.5 were used (most of the energy dissipated in the active corona). The contours represent the observed correlations 
(Zdziarski et al. 1999; Gilfanov et al. 2000). Panel b) shows hard X— ray flux - spectral index anti-correlation. Each "triplet" of points 
corresponds to a given value of rj and f = 0.1, 0.3, 0.5. The anti-correlation, which has an obvious physical interpretation in this model, 
seems to be opposite to what is observed in the data (e.g. Done et al. 2000). 



case R{C) ~ 0.5 when the anti-correlation between R and 
C{t) can hold through the entire range of C{t). Details of the 
relation 7?(C)-Aavar depend on the mapping C{t) rf), 
i.e. which parameter, ^, rj or both of them, cause flux vari- 
ations. Solid squares in Fig. |^ show the case of variations 
C{t) caused by changes of ^ only (for 77 = 1). Reduction 
of the line variability is strongest in this case, since a given 
variation of C{t) — Fx cau ses l argest possible change of R, 
as already discussed in Sec. 3.3 and shown in Fig. 0. On the 



other hand, the available range of R(C) is rather limited. 

These results are quantitatively dependent on the as- 
sumed constancy of X-ray source geometry (see comments 
following Eq. (M)). The reduction of variability, AcTvar, 
could be larger, if the geometry were correlated with the 
X-ray luminosity, changing in such a way as to make the 
illuminating flux. Fx, to increase more strongly than the 
observed flux, Jhard- However, as we are going to argue in 
the next Section, the combined timing and spectral obser- 
vational result put quite tight constraints on any geometry- 
luminosity dependences. 

Finally, we note that both R{C) and (7var(F) — o"var(i^ 
are observable quantities and their relation plotted in Fig, 
can be used to test the model. 



5 SPECTRAL DATA 

5.1 The J?— F— flux correlations 

The spectral index of the hard X-ray power law, F, and the 
observed amplitude of the reprocessed component, R, ap- 
pear to be correlated in the data of both AGN and Galactic 
Black Holes (GBH) (Zdziarski, Lubihski & Smith 1999; Gil- 
fanov, Churazov & Revnivtsev 2000): the harder the spec- 
trum (smaller F), the smaller the amplitude. Two mod- 
els were considered to explain the correlation: truncated 
cold disc with hot central flow (e.g. Esin, McClintock & 



Narayan 1997) and relativistic outflow of emitting plasma 
(Beloborodov 1999ab). The model of an accretion disc with 
a hot layer is also able to qualitatively explain the correla- 
tion, with Thot being the control parameter: the larger rhot 
is, the smaller the effective R and the harder the spectrum, 
due to reduced soft flux from thermalization (assuming that 
most of the energy dissipation occurs in the active corona) . 
We present here the -R-F correlation predicted by the model. 

(equation 



The estimate of R is obtained as in Sec. 3.1 



while for F we use the formula from Beloborodov (1999ab), 



2.33(/h//s - 1) 



-1/10 



(15) 



where /h and Is are the power in hard and soft radiation 
crossing the active region, respectively. We estimate 



(16) 



= y [Fdissip(At) + -Frcproc(At)] d^i 

where Fa^sip = Fu = C^tot and 

Froproc = (l-a)Jk = i?7(l-C)(l-a)-FtotPtr (rhot(C,'?))-(17) 

The second term describes the contribution to the soft flux 
from the thermalized fraction of the illuminating X-rays. 
We have also introduced the geometrical factor /j.^ (see Be- 
loborodov 1999ab), describing the fraction of soft flux ac- 
tually intercepted by a magnetic flare (in computations we 
assume isotropic emission). The hard power is simply esti- 
mated as /h = 2Fx = r?(l - C)-Ftot. 

Figure |^ presents results of the computations for a num- 
ber of /is values. The shape of the model correlation is some- 
what different from the correlation observed in the data. 
The model requires /is > 0.6 in order to obtain quantitative 
agreement with the data, i.e. 1 — fia < 0.4 of the soft flux 
crossing an active region. It has to be remembered however 
that the observed correlation, particularly among different 
objects, is more likely to be due to changes of geometry 
driven by e.g. m rather than varying output of magnetic 
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activity, and thus parameterizing the correlation by our pa- 
rameters ^ and r) may not be sufficient. 

Another correlation applicable to AGN observations 
may be the -Fhard-r (anti)correlation plotted in Figure It 
follows from combining results shown in Fig. ^ with those in 
Fig. 0: an increase of illuminating flux, Fx, gives a decrease 
of R and a decrease of F. If the observed X-ray flux, -fhard, 
follows the illuminating flux, fx, (i.e. if the Eq. (^) holds), 
the predicted -Fkard-r anticorrelation is obtained. However, 
an opposite effect seems to be observed in the data from 
RXTE observing campaigns of some AGN, where the spec- 
trum generally softens (F increases) when the observed X- 
ray flux goes up (NGC 5548, Chiang et al. 2000; IC 4329a, 
Done et al. 2000). At least in IC 4329a the spectral analysis 
of Done et al. (2000) reveals that this is not an artifact of lim- 
ited band pass of RXTE instruments. If this observational 
result is generic, it would mean changes of the flares' ge- 
ometry correlated with their luminosity, so that the relation 
between Fhard and -Fx can be reversed. From the discussion 
of variability (Sec 4.S), we need a correlated increase of rhot 
with the flare luminosity. If this is to be accompanied by a 
softening of the spectrum, the flare's geometry has to change 
in such a way, that the intercepted soft luminosity is larger 
than the increased hard luminosity (Zh/^s has to decrease; 
Eq. (p^), despite the decreased soft flux from thermalized 
hard X-rays. 



5.2 Model dependence of relativistic smearing 
and amplitude of reflection 

The hot layer is thickest at close distances from the central 
black hole, accor ding to results of computations of rhot('") 
presented in Sec. 3.2. The line generation is then sup- 
pressed at small r, which means that the inner disc radius de- 
termined from the line proflle may be larger than rms = 6 rg. 
We quantified this effect in our model, since indeed the line 
observed in a number of sources is narrower than in the ex- 
treme case of MCG-6-30-15 (IC 4329a, Done et al. 2000; 
Cyg X-1, Done & Zycki 1999; GS 2023-H338, GS 1124-68, 
Zycki, Done & Smith 1997, 1998). 

We computed the relativistic smearing in the repro- 
cessed spectra for the radial dependencies Thot('") plotted 
in Fig. M. We assumed that the irradiation emissivity fol- 
lows radial distribution of gravitational energy dissipation, 
G(r) (equation Ratios of the reprocessed spectra to the 
reference case of no hot skin (so full relativistic smearing), 
are plotted in Fig. ^. The magnitude of the effect obviously 
depends on the thickness of the HL. Even for the thick- 
est HL considered (^ = 0.01, maximum rhot('") ~ 0.5), the 
"wiggles" in the ratio plots are only a few percent, for pure 
reprocessed component, i.e. no contribution from primary 
photons. It would then be rather difficult to detect such an 
effect in the data, the more so that the additional broaden- 
ing due to Comptonization in the HL (not included here) , is 
strongest for the thickest HL. Qualitatively then, reduction 
of effective amplitude of the reprocessed component is much 
more prominent, than any changes to the proflle of the Kq, 
line (as is also clear from Fig. ^). 

One can estimate the maximum magnitude of the sup- 
pression of the broadest line component, noticing that it will 
occur for rhot('") being a step function; Thot very large for 
radii r < r^ut, and rhot ~ otherwise. In other words, no 
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Figure 8. Reprocessed spectra computed for the radial Thot('') 
dependencies plotted in Fig. ^ (solid lines), compared to the ref- 
erence case of no hot skin (full relativistic smearing; dashed line). 
Since Thot('') reaches maximum at r ^ lOrg, the most strongly 
relativistically smeared component is suppressed in the spectra. 
Upper panels show that the main effect is a decrease of the overall 
amplitude, described by Eq. Lower panel show the difference, 
ionized skin spectrum — reference spectrum, (after correcting 
for the decrease of amplitude), normalized to the primary power 
law spectrum. The additional Comptonization broadening of the 
spectral features is not included. 



observable reprocessed component is produced for r < r-cut 
(Ross, Fabian & Young 1999) . This obviously means a reduc- 
tion of the effective amplitude of the reprocessed component, 
which one can relate to Tcut 

To do so, we write R — i^rcfl/(aio-fdiroct), i.e. R is the 
reflection albedo (ratio of flux of reflected radiation to the 
flux directed towards an observer) , normalized to the albedo 
ao of cold reflection, expected from an isotropic source above 
a flat disc. Both fluxes can be written as radial integrals, e.g. 



Frofl = 27r 



Frca{r)rdr. 



(18) 



Averaged over many dynamical time scales, the radial dis- 
tribution of energy generation in the active corona can be 
expected to follow the gravitational energy dissipation pre- 
scription. If isotropy of emission is assumed, then 



i^rcfl (0 =-Fdircct(r) =G(r). 

If a purely scattering layer of a thickness giving ; 
sion probability pti (reflection probability 1 —pti] 
reflecting disc up to a radius rcut, then 

Frcti{r)/ao = 27r plG{r)r dr + 

2tt G{r)r dr 

and 



(19) 

I transmis- 
covers the 



(20) 



fdircct(r) = 2n J^^2 (2 - ptr)G{r)rdr + 

2tv /J'""" G{r)rdr. (21) 

The resulting function R{rcut) is plotted in Figure ^ for 
Ptr ~ 0, which maximizes rcut for any R. For ptr('") — 1 
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Figure 9. Relation between effective amplitude of the repro- 
cessed component and the amount of relativistic smearing, for 
an ad hoc model, where the disc is assumed to be covered by a 
perfectly reflecting HL (i.e. Tj^ot S> 1) from rms up to a radius 
Tcut- This reduces the effective amplitude of reflection and re- 
moves the broadest component from the spectral features, with 
?'cut being the inner disc radius as inferred from fitting the profile 
of the Ka line. 



(no hot layer) we obviously obtain i?(rcut) = 1, but for any 
Ptr(f) < 1 the amplitude R decreases dramatically with in- 
creasing Tcut, i.e. when the ionized zone becomes more ex- 
tended. The reason for the dramatic decrease is that the 
emission is strongly concentrated towards the center. Inter- 
estingly, many Galactic black hole binaries show i? ~ 0.3 
and Rin ~ 20r-g, in accord with the relation shown in Fig- 
ure ^ but we note that, with the assumption of perfect re- 
flectivity, the geometry here is equivalent to the disc being 
truncated at rcut , and replaced by central X-ray source. We 
also emphasize that the assumption pti{r < rcut) = 0, is not 
justified by our present understanding of the model. 



5.3 The X-ray "Baldwin effect" 

The observed equivalent width (EW) of the Fe line ap- 
pears to be anticorrelated with average 2-10 keV luminosi- 
ties of the sources (Iwasawa & Taniguchi 1993; Nandra et 
al. 1997), in a similar manner as EW of CIV A1550 line is 
anticorrelated with the UV continuum luminosity (Baldwin 
1977). This could be explained by the presence of the hot 
layer, since indeed Thot is larger for larger Fx, although in 
this case we expect the flux to be driven by accretion rate, 
rather than ^ and/or rj. Since the average X-ray luminos- 
ity can be expected to be proportional to the accretion rate 
M — rhMEdd oc rhM, we have computed the EW varying 
either M ox m and keeping the other parameter constant. 
We again assumed that situation at r = 10 Tg is character- 
istic to the whole disc, and we adopted ^ = 0.5 and 77 = 1. 
Results presented in Fig. |l^ show that indeed the EW of the 
Kq line (oc i?) decreases with m. The reason for the under- 
lying increase of Thot with rh is that the increase of with 
rh is slower than linear. The gravity parameter A oc Hg^/m 
then decreases with rh. We note that the vertically averaged 
SS solution predicts i/d oc rh and, in consequence, constant 
Thot in the regime of Pj-ad S> -Pgas and k ~ Kcs. In our ver- 
tically explicit solution Pgas is not neglected. Contribution 



M=10' M, 
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Figure 10. Equivalent width of the Fe Ka line (solid curve, left 
scale) as a function of the average X— ray source luminosity, Lx, 
assuming that Lx oc m, i.e. M = const. EW decreases with m 
which could explain the observed X— ray "Baldwin effect". The 
dotted curve shows the corresponding gravity parameter A, which 
decreases due to slower than linear increase of the disc thickness 
with rh. Parameters: M = 10^ Mq, « = 0.1, § = 0.5, 77 = 1. 

from convective energy transfer contributes to departures 
from the SS solution as well. 



6 DISCUSSION 

Presence of a hot, ionized layer on top of an X-ray illu- 
minated accretion disc has inevitable observational conse- 
quences. Most importantly, it reduces the amplitude of the 
effective, observed "cold" reprocessed component. Since the 
thickness of the hot layer (hence the reduced amplitude) is 
correlated with the X-ray flux, a number of observable cor- 
relations follow from this model. 

Generally, the predicted Fe fluorescent line originat- 
ing from the disc is weaker (relative to the continuum) when 
the source is brighter (provided that the illuminating spec- 
trum has r < 2 so that there is no He-like Fe in the hot layer; 
NKK) . This effect can be used to explain two observed phe- 
nomena. Firstly, the reduction of the amplitude of variabil- 
ity of the line with respect to the driving continuum {E > 7 
keV) in AGN monitoring campaigns (Reynolds 2000; Done 
et al. 2000; Chiang et al. 2000). Quantitatively, our com- 
putations indicate that it is not possible to obtain constant 
absolute line flux constant. That is, the relative amplitude of 
the reprocessed component decreases more slowly than Px ^• 
Nevertheless, reduction of the variablity amplitude by a few 
per cent is possible. However, the variability reduction is in- 
versly related with the mean amplitude of the reprocessed 
component (i.e. observed at the mean flux level). Thus, in 
sources like MCG-6-30-15, where P ~ 1, this effect should 
actually be negligible. Moreover, the X-ray continuum in 
MCG-6-30-15 is rather soft (F ^ 2), which implies a contri- 
bution from a highly ionized reprocessed component (NKK), 
that does not seem to have been observed in this source. 
Perhaps the hot skin developes on accretion discs in the 
other Seyfert 1 galaxies, where both P < 1 and ^constant 
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Fe Ka line flux was observed (IC 4329a, Done et al. 2000; 
NGC 5548, Chiang et al. 2000). Quality of current RXTE 
data should be good enough to quantitatively test the model 
relation presented in Fig. H. 

Second, the decrease of the EW of the Fe Kc, line in 
brighter sources (X-ray "Baldwin effect" ) can be explained, 
if the inverse relation between efficiency of the line produc- 
tion and X-ray flux is applied to a sample of objects. This 
is not quite a trivial statement, since different average X- 
ray luminosities of sources are rather due to different mass 
accretion rates, than a variable output of some magnetic 
activity. 

The model is also able to explain the correlation be- 
tween the amplitude of the reprocessed component and X- 
ray spectral index, again observed both in a sample of ob- 
jects (Zdziarski et al. 1999), and in the time history of a 
single object (GX 339-4, Gilfanov et al. 2000; GS 1124- 
68, Zycki et al. 1998). Our computations of the correla- 
tion give perhaps somewhat different shape to what is ob- 
served, although the observed correlation is also subject to 
instrument -related uncertainties. 

Another relation often looked for in data and used to 
constrain theoretical models is a relation between spectral 
index and hard X-ray flux. This is usually more difficult to 
obtain from data because of uncertain corrections from lim- 
ited detector band pass to the entire X-ray band. Perhaps 
the most robust result was obtained by Done et al. (2000) 
for IC 4329a, where broad band spectral analysis revealed 
that the spectrum became harder when the X-rays {E > 1 
keV) decreased. Similar relation was reported for NGC 5548 
(Magdziarz et al. 1998; Chiang et al. 2000) aflhough this in- 
volved flux in 2-10 keV band. On the other hand the soft X- 
ray transient GS 1124-68 (Nova Muscae 1991) showed a clear 
hardening of its spectrum during its decline phase (Zycki et 
al. 1998; Zycki 2001), which was however correlated with an 
increase of the E > 1 keV flux by « 50 per cent (so called 
're-flare'). The model discussed in this work predicts anti- 
correlation of r and -fhard, i.e. hardening of the spectrum 
when a source brightens, unless there are correlated changes 
of the X-ray source geometry and its luminosity. 

The structure of the hot layer and resulting correla- 
tions between observables may also depend upon the model 
adopted for the underlying accretion disc. In this paper we 
adopted the aPtot SS prescription for viscosity. However 
such discs are subject to a number of instabilities in the pa- 
rameter range appropriate to AGN. Classical theory predics 
thermal and viscous instability whenever contribution of gas 
pressure to Ptot drops below 0.4 (Kato, Fukue & Mineshige 
1998), which in practice means the region within a few hun- 
dreds rg, for rh a few per cent. Therefore, models with a 
different viscosity prescription were sometimes considered, 
e.g. a-Pgas, which are free from those instabilities. Solving the 
vertical structure of aPgas discs for a number of (5, rj) values 
we flnd that the disc height is very similar (less than 1 per 
cent difference) to that of the corresponding aPtot disc. The 
resulting relations between observable quantities are then 
practically identical to those for aPtot discs. Our current 
understanding of mechanisms providing viscosity strongly 
points out towards the magneto-rotational instability as the 
source of turbulence responsible for the viscosity (see Balbus 
& Hawley 1998 for review) . Global magnetohydrodynamical 
simulations reveal the disc structure very different from the 



smooth, classical, hydrodynamical models (Hawley 2000). 
The spatial and temporal structure is highly turbulent and 
inhomogeneous, displaying fractal-like scalings (Kawaguchi 
et al. 2000). The relevance of the presented model to such 
inhomogeneous, clumpy discs needs further studies. 

The presented model is currently one of at least three 
general scenarios within which it is possible to explain bulk 
of observed data from low/hard states of AGN and GBH: the 
other two are the cold disc with inner hot flow model (Esin et 
al. 1997 and references therein; Rozariska & Czerny 2000b) 
and relativistically outflowing plasma model (Beloborodov 
1999ab). None of them fully addresses all complexities of 
observed spectral and timing data, but all three account for 
the main facts: the range of observed spectral indices, the 
R-F correlation, relativistic smearing of Fe reprocessed fea- 
tures (see Poutanen 1998 for review and e.g. Di Salvo et 
al. 2001 for application to Cyg X-1). It also seems possi- 
ble to incorporate into all of them the basic X-ray timing 
observables: the shape of PSD, time lags, auto- and cross- 
correlation functions (Poutanen & Fabian 1999; Bottcher 
& Liang 1999; see review in Poutanen 2000). Clearly, fur- 
ther progress can only be made if the models are developed 
in more details incorporating more relevant physical pro- 
cesses. Confronting detailed, quantitative predictions with 
good quality spectral and timing data should hopefully en- 
able distinguishing which, if any, of our current ideas are 
closest to reality. 



7 CONCLUSIONS 

• The model does predict certain reduction of the vari- 
ability of the Fe line, compared to the variability of its 
driving continuum, although it does not seem possible to 
obtain absolutely constant line flux. 

• The anti-correlation of the equivalent width of the Kq 
line and average source luminosity (the X-ray "Baldwin ef- 
fect") is reproduced in the model. 

• Reduction of the amplitude of the observed, "cold" re- 
processed component is more evident than any reduction 
of the spectral broadening (due to relativistic effects and 
Comptonization) of the Ka line. 

• The model is able to quantitatively explain the observed 
R-F correlation. 

• The geometry of the X-ray sources is quite tightly con- 
strained, if the model is to explain both timing and spectral 
observational data. 
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